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ION CYCLOTRON WAVE GENERATION IN UNIFORM AND NONUNIFORM 
PLASMA INCLUDING ELECTRON INERTIA EFFECTS 
by Donald R. Sigman and John J. Reinmann 
Lewis Research Center 

SUMMARY 

Theoretical calculations were made for radiofrequency power transfer to ion cyclo- 
tron waves from a Stix coil. The analysis includes electron inertia terms and is applied 
to a cylindrical column of cold, collisionless atomic- hydrogen plasma with and without 
radial density variations. The results showed both an optimum coil wavelength and an 

optimum density for efficient power transfer. For typical plasmas, the optimum density 
12 12 

is between 1x10 and 5x10 per cubic centimenter. Without a radial density gradient, 
the inclusion of the electron inertia terms reduced the predicted power transfer. Several 
radial density distributions were considered. Power transfer to inhomogeneous piasmas 
is maximum when the volume average density is near the optimum density computed for 
homogeneous plasmas. Curves of power absorption as a function of magnetic field for 
inhomogeneous plasmas were similar to those for constant density. Under some condi- 
tions, radiofrequency power is coupled predominantly to a single plasma mode of oscilla- 
tion; under other conditions, the power may be distributed among several modes. Thus, 
measurements of plasma wave properties may prove difficult to interpret. 


INTRODUCTION 

The thermalization of ion cyclotron waves in a plasma column offers a promising 
technique for heating plasma ions to extremely high temperatures. Stix (ref. 1) devised 
a generating structure (radiofrequency coil) to couple radiofrequency power to ion- 
cyclotron waves efficiently and has also made a careful theoretical analysis of the power- 
coupling process. In his theory, Stix assumed a cold, collisionless magnetoplasma and 
neglected electron inertia terms in the equations of motion. In reference 2, Stix’s theory 
was employed to determine the effect of various plasma and radiofrequency coil param- 
eters on the power transfer to ion cyclotron waves in a homogeneous, atomic-hydrogen 


plasma. In the present report, the parametric analysis is repeated, but electron mass 
is not assumed to be zero. Furthermore, both homogeneous and radially inhomogeneous 
plasmas are considered. The results of this report were first presented at the 1966 
Summer Meeting of the American Physical Society, Minneapolis, Minnesota. 

As in references 1 and 2, the radiofrequency coil is represented by an azimuthal cur- 
rent sheet of finite length wrapped around an infinitely long plasma column that is placed 
in a vacuum. The plasma is immersed in a steady and uniform magnetic field parallel to 
the axis of the column. 

In the region underneath the finite- length current sheet, power is transferred to the 
ion cyclotron wave. This wave propagates unattenuated out both ends of the cylindrical 
current sheet and down the plasma column. Since the plasma model is cold and collision- 
less, all energy from the current sheet is converted to electromagnetic energy in the 
wave. 

The theoretical analysis proceeds as follows: Maxwell's equations are employed with 
the equations of motion for the plasma ions and electrons in order to derive the fourth- 
order partial differential equations for the electric and magnetic fields in the plasma. 
Then, expressions are obtained for the fields under infinitely long current sheets with 
wave number k. Following a Fourier integral technique used by Stix (ref. 1), one uses 
the fields derived for infinitely long current sheets to obtain the fields produced by a 
finite-length current sheet. For the homogeneous plasma, a dispersion relation is de- 
rived. 

The equations derived herein were used to make digital computer calculations to de- 
termine the effects of coil wavelength, overall coil length and radius, magnetic field, 
plasma density, and plasma radius on the power transfer to ion cyclotron waves in a 
homogeneous plasma. In addition, sufficient calculations were made for inhomogeneous 
plasmas and for an annular cylinder of plasma of uniform density to determine the major 
effect of radial density variations on power transfer. 

The assumption of zero electron mass simplifies the problem considerably. Early 
attempts by the authors of this report to extend the analysis of reference 2 to inhomogene- 
ous plasmas resulted in the prediction of unrealistically low power-transfer efficiencies 
for slight departures from uniformity. Furthermore, the findings of reference 2 yield 
the unrealistic result that maximum power transfer monotonically increased as the den- 
sity decreased. The inclusion of electron mass in the analysis removed both these dif- 
ficulties. 


THEORY 

To derive equations for radiofrequency power transfer to a plasma, a geometric con- 
figuration similar to that used by Stix (ref. 1) is assumed. The elements of the config- 
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Vacuum 



''-Azimuthal current-sheet density, j*e'* kz " 

Figure 1. - Plasma column and current sheet. 


uration are shown in figure 1. A cylinder (solid or annular) of cold plasma, infinitely 

long, is assumed to be surrounded by vacuum. A uniform, static magnetic field B Q 

(maintained by external coils) is parallel to the z-axis of the cylinder. The outer radius 

of the plasma is p, and the inner radius is £. For a solid cylinder of plasma, the inner 

radius £. is zero. At the radius r = s (where s > p), there is an azimuthal current 
i(k z-wt) 

sheet of density j*e ° , where k Q is the current- sheet wave number for the z- 

direction. Because the current sheet for the case of interest is of finite length, it is nec- 
essary to apply a Fourier integral analysis (ref. 1) to express the finite- length sheet as 
an integral over k of a continuous spectrum of azimuthal current- sheet components, 
each of which can be considered to be of infinite length. In this approach, the equations 
are formulated for an infinite -length current sheet and then are combined in the manner 
indicated by the Fourier integral analysis to obtain the results for the finite sheet. 


Equations of Motion and Maxwell's Equations 

For an infinitely long current sheet, the time-varying electric and magnetic fields 
and the currents and particle velocities are Fourier analyzed in space and time according 
to and will have the general form 

F(r, z, t) = f(r)e‘! kz ' a,t l (1) 

All terms are assumed to have azimuthal symmetry. For a cold, collisionless plasma, 
the appropriate equation of motion for ions is 


m ; 



+ 


li 



v, x B 

= eE + e 

c 


(2a) 


and for electrons, 


(2b) 


9v v X B 

m — - + v - Vv = -eE - e 

e 3t “ e _e ~ c 

(All symbols are defined in appendix A and cgs-Gaussian units are used throughout this 
report. ) A Fourier analysis is made in time according to equation (1). No zero-order ve- 
locities are assumed and B = B Q + B(r). Then the higher order v. Vv and v x B(r) 
terms are dropped, and equations (2a) and (2b) become 


[v, (r) X B 1 

-iwmjV^r) - eE(r) + e L J 

c 

(2c) 

[v (r) x B J 

-iwm e v e (r) - -eE(r) - e ^ e bJ 

(2d) 


c 


The plasma is assumed electrically neutral: n 0 = n^ = n(r). The current density in the 
plasma is given in terms of the macroscopic particle velocities as 

_j(r) = n(r)e|y.(r) - v g (r)J (3) 

Equations (2c), (2d), and (3) are used to eliminate v. and v g to obtain _j(r) as a function 
of E(r) 


_j(r) = — — K’(r) -E(r) (4) 

47ri 

The tensor K'(r) is given as 



in the coordinate system where the electric field is given as 
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R’ = 


L’ = 




S' = - (R' + L’) 
2 


D’ = - (R' - L’) 
2 


(5c) 

(5d) 

(5e) 

(5f) 

(5g) 


where Il e and IIj are the electron and ion plasma frequencies, and and are 
the electron and ion cyclotron frequencies. Maxwell’s equations are written as 


« 9B 

V x E = — (6a) 

_ c 9t 


4 Jr 1 OCj 

V X B = — _j + - — (6b) 

c c 9t 


and are combined to give 


V X V X E = 




(6c) 


where the Fourier analysis is again taken in time. Now, substituting equation (4) into 
equation (6c) yields 


V X V X E = 


— K(r) • E 
c 2 ~ “ 


(7a) 
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where K = K' + 1. The three components of the vector equation (eq. (7a)) are 


’3E \ 2 

ik | — - - ikE ] - — (SE - iDE J = 0 
3r 7 c 2 


(7b) 


k 2 E fl rEJ - (iDE + SE fl ) = 0 


dr \r 3r 7 2 


(7c) 


i _a_ 
r 3r 


r I ikE 


3r 


_ PE = 0 

2 z 
c 


(7d) 


where the Fourier analysis is now taken in the z-direction; E r , E^, and E z are func- 
tions of the radius r only; and S = S’ + 1, D = D', P = P* + 1, R = R f + 1, L = L' + 1. 


Boundary Conditions 

The boundary conditions that are to be imposed on the solutions to Maxwell’s equa- 
tions are as follows: 

Plasma- vacuum- interface: The components of the tangential E-fields and B-fields 
must be continuous across a plasma-vacuum interface (r = p and r = £): 


or 


^ E plasma\ 

^ E plasmaj 


1 __ 

(El aC 

| 

r=p,£ 

V e ) 

r=p,£ 

r=p,£ 

/■nV ac N 

5 > 

L=p,£ 


(b™ c ) 

1 

r=p,£ 

\ 6 ) 

1 

II 


ikE 


aT? plasma\ 
plasma dhj z 1 


= ikE 


3r 


’ r=p, £ 


3E vac ' 
vac z_ 

3r 




' r= P> £ 

j 
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B; 


plasma\ 


Jr, 


r=P ,JL 



£ 


or 




'3E? laSma i 
d 

dr 


r 9E 


vac\ 


f r=p, £ 


3r 


' r =P» !■ 

J 


Current sheet at r = s: The components of the tangential E- fields and Bg must be 

continuous at/ r = s. ,The component B_ is discontinuous by an amount (4jr/c)j*: 

z 


/ B plasma\ _ / B vac\ _ j 
\ z / r=s V z / r=s c 


or 




r 3E^ lasma i 

dr 


r 3E 


vac\ 


'r=s 


3r 


47riw . * 
3* 


r=s 


Vacuum Solutions 


The solutions to Maxwell’s equations in the vacuum regions (with displacement cur- 
rent neglected) are 


Eg = bjIjCkr) 

E z = Cl I 0 (kr) 

E r = -iCjIjCkr) 

* C l I l (kr) 

c kc 

i - B z = kbjBCkr) 
c 


(8a) 

(8b) 

(8c) 

(8d) 

(8e) 
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E e = djljfkr) + dgKjfkr) 
E z = ei K 0 (kr) 


E r = 


pSr^ s 


iZBg^eftQ*) 

c kc 2 


i - B z = kdjI/jCkr) - kd 2 K Q (kr) 

c J 


E 0 = fjKjCkr) 


E z = * 1 * 0 <& 
E r = iejKjOkr) 


i- B e =^ e l K l^) 


i - B z = -kf^kr) 
c 


where bj, c«, dj, d 2 , e^, and f ^ are constants of integration. 

For the solid plasma cylinder (£ = 0), the solutions for Olr^t are ignored. 


Homogeneous Plasma (Solid or Annular) 


For constant density (n(r) = Const), equations (7b) to (7d) are combined to give 


(ST - + C)E fl = 0 
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(9b) 



ikPE z = - 


1 _d_ 

r 9r 


r(SE r - iDE 0 ) 


3/19 \ 

where if is an operator given by — [ r ] 

3r \r 3r / 


k 2 (P + S) - — (PS + RL) 
B _ k 2 (P + S) co 2 p | a) 2 (D 2 - S 2 ) c* 

S c 2 c 2 S S 


k 4 - 2 — k 2 S + — (S 2 - D 2 ) 

_ p 

? k 2 - < 2 ! R Vk 2 - 4 

2 4 

s 


c C 


\ c / \ c / 


(9c) 


For constant- density plasmas, B and C are not functions of the position variable r, 
and the fourth-order differential equation (9a) is factored into two second-order Bessel 
differential equations 


if 


2 

+ V 1 


E„ = 0 


(10a) 


if + v‘ 


E e = 


where 


v 


2 

1 



v 


2 

2 “ 



4C 


(10b) 


(10c) 


(lOd) 


From equations (10a) and (10b) four Bessel function solutions are obtained for E n : 

C7 

Jl(^l r ), Yj(^jr), Jj^r), and The complete solution for E^ is then 
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E e = a l J l^l r ) + a 3 Y l^l r ) +a 2 J 1 ^ 2 r ) +a 4 Y l^2 r ^ 


for $_ S r ^ p 


(11a) 


2 2 

where a^, & 2 , a^, a^ are constants of integration. (It is assumed that v ^ and are 
both positive. When either v 1 2 or r 2 becomes negative, the appropriate Bessel func- 
tions should be substituted. See also appendix B.) From equation (9b) 


iD— E = fk 2 - SC- — SlE 


e 


n Z' 2 or c , 
= (k + - — S] 

1 2 
c 


+ agYj^jr) 


+ 




a 2 J 1^2 r) +a 4 Y l 



is obtained, and from equation (9c), 


ikPE 

z 


i _d_ 
r 3r 


r(SE r - iDE 0 ) 


and 



av+ 


a l J 0 (^i r ) + a 3 Y o( y l r ) + & v 2 a 2 J 0^2 r ^ + a 4 Y 0^2 r ^ 


1 

are obtained. From V x E 

~ c 3t 


(lib) 


1 c Bz = r ~t ^ = y l[ a l J 0 (i; l r) +a 3 Y 0^1 r ) 


+ V t 


a 2 J 0^2 r ) + a 4 Y 0^2 r ^ ( 12a ) 


and 
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i — B„ = ikE, 


aE„ 


-= a 


3r 


a l J l(^l r ) + a 3 Y l^l r )] + ^[ a 2 J 1^2 r ^ + a 4 Y 1^2 r ^] ( 12b ) 


are obtained. Equations (Ha), (lib), (12a), and (12b) are for f ^ r ^ p, where 


o i, 2 2\ DI id 

S(k + r ^ 1 - RL — 


a = 


kD 


S (k 2 + v\ 


- RL 


c o 


P = 


kD 


The constants a-. , a 0 , a„, a. must be determined from boundary conditions. For the 

^ plasma 

solid cylinder of plasma, 1 = 0, and a^ and a^ = 0 because E^ is finite at r = 0. 


Power Transfer to Homogeneous Plasma 


When the aforementioned boundary conditions are applied to the constant- density solu- 
tions (solid or annular), a system of simultaneous equations is obtained for the coeffi- 
cients in equations (8), (11), and (12). These simultaneous equations can be used to solve 
for ap ag, a^, a^, dp d ep and fp Calculation of the power absorbed by the plasma 
requires the value of the azimuthal electric field at r = s (ref. 1). This value can be 
written as 


En aC (k) 


r=s 


= f^^ks) 


(13) 


The electric field 


E^ aC (k) 


is the field at the position of the current sheet for an in- 


r=s 


finitely long current sheet. What is actually needed is an expression for (E^) for a 


finite- length current sheet of length 2a with a current density 

i(k 0 z-u>t) 


r=s 


j* = j*e 


(-a iz? a) 


(14a) 
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y * = 0 | z | > a 


(14b) 


For a finite- length current sheet, (E„) is obtained as follows (ref. 1): The 

y r=s 

current- sheet density _j* from equations (14) is decomposed into its Fourier integral 
spectrum, so that 


j*= ej* 


/ 


Slp(k - k o )a e i0=— t) ^ 


7r(k - k Q ) 


(15) 


The electric field for the finite-length current sheet can then be represented as the sum 
(integral) of the electric fields of a continuous set of infinitely long current sheets, each 
having the infinitesimal amplitude 


_ .tpa-Mt) 

7T(k - k o ) 


dk 


so that 


(E„) 


r=s 




sin(k - k Q )a 
7r(k - k Q ) 


e i(kz-wt) 


dk 


(16) 


Equation (16) was integrated by contour integration (as in ref. 1) to give the following re- 
sult for a right-running current- sheet wave: 


(E fl ) = 3* 


r=s 



ik Q z-iu)t 


n=l 


1 - e 


- i ( k o" k n )( z+a ) 


k o" k n 


-a< z < a 


(17a) 


where 
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(17b) 


R n 


(k-k n ) 



r ^k=k. 


n 


j* 


(It should be noted that the R n defined in eq. (17b) is not the same as the R n defined on 
p. 96 of ref. 1. ) 

The values k = k are the values of k where 
n 

correspond to the poles of the integrand in equation (16). The k n 's are the wave num- 
bers for the natural modes of the plasma column in a vacuum. This suggests that each 
term in the sum in equation (17a) is the electric field of the n 1 natural mode. 

The component of the electric field (E fl ) , which stems from power absorption in 

w r=s 


nvac\ 

J0 / r=s 


becomes infinite. They 


the plasma, tends to decelerate charge carriers in the current sheet. Thus, the instan- 
taneous power per unit area required to drive the surface current against the electric 


field at r = s is 


[V>*] 


Therefore, the time-averaged power transfer to the 


r=s 


plasma is obtained by integrating over the total surface of the current sheet as follows: 


P t = - Re 
r 2 


£ 


- (E fl ) 


-ik z+icot 


j*e 


2tts dz 


r=s 


(18) 


Integrating equation (18) gives 


or 


P t = (27ras)(-iaj* 2 ) ^ R n S r ^n^ (!9a) 

n=l 


QO 

P t = (27ra 2 )j* 2 ^ R' n S r (v n ) (19b) 

n=l 


where 


2(1 - cos 7] ) 
S^ n ) = 7 ~ 


V 


n 


(19c) 
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= 2 < k o - k n> a 


(19d) 


and 


R' = -isR 
n n 


(19e) 


From equations (13) and (17b) and from the fact that the denominator of 
(hereinafter denoted as Tj(k)) must be zero, R^ can be calculated. 



r=s 


R' n for Homogeneous Plasma 

For the simpler case of the solid cylinder (SL = 0) the expression for R^ is 



s 2 K^(ks) 


k P K Q (kp)K 1 (kp) 


|a[pJ 1 (i/ 1 p) + yy 1 pJ 0 (i/ 1 p)]j 1 (p 2 p) - /3[pJ 1 (i^ 2 p) + yy 2 pJ 0 (i^ 2 p)]j 1 (i/ 1 p)J 


rj(k) 


( 20 ) 


k =k n 

r i=° 


where 


Kj(kp) 

kpKg(kp) 

The functions r^(k) are given in appendix B. The indeterminancy in equation (20) may be 
eliminated by L' Hospital's rule, so that 


k “ K 

lim 5 — = 

k-k r l (k) 

K K n 1 k=k 

a 


dr,(k) 


dk 


Jk=k 


n 


( 21 ) 


For the case of the hollow cylinder of plasma, 


4jho\ s 2 [k 1 (ks)] {[~ M 41 J 1^1 p ^ + M 42 J 1^2 p ^ ' M 43 Y 1^1 P ^ + M 44 Y 1^2 P ^]} 


R; = (k * k J 


2 J kpK 0 (kp)K 1 (kp) 


r 9 (k) 


k=k n 

r 2=0 


( 22 ) 
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I 


where ^(k) = det M; and M^, M^g, are the minors of the matrix M, which 

is also given in appendix B. Again, the determinant M becomes zero when k = k n and 
the indeterminancy is eliminated by L’ Hospital's rule. The values v^, and k R , 

when Tg = 0, are the values obtained for the natural modes of oscillation of an infinitely 
long plasma column surrounded by a vacuum. Thus each term in the sum of equation (19b) 
may be throught of as the power transmitted by a single natural mode of the plasma. 


Inhomogeneous Plasma 

If, as in an inhomogeneous plasma, n(r) does not equal a constant, equations (7b) to 
(7d) combine to give a lengthy, unfavorable fourth- order differential equation for E^, 
which must be solved by numerical means. If the 1 is neglected in the terms R, L, S, 
and P (which is the same as neglecting displacement current), R, L, S, P, or D can 
be written in the form Q = Q Q g(r), where Q q is the value of R, L, S, P, D, or n at 
r = 0, and g(r) is the density distribution function so that n(r) = n Q g(r). Thus, when dis- 
placement current is neglected, the differential equations for the electric- field compo- 
nents in an inhomogeneous plasma are 


d \ , d \ fZ g'\ , 3 \ 

dr 4 dr 3 \ r «/ 3r 2 


_ 3 _ _ 2 g[ + P _o 

2 r g S 

r & o 



-k 2 S + — 

O r 


"O L_ 


S^-D 2 )g 


3E, 




3r 




r 








+ V 

V _o 

s 



0 

V- 

\ c / c 



+ ik 2 £ + i? 

r g 2 
& r 


"A 



J 


(23a) 
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(23b) 



dr 






\T" 






J 


E 


z 


-c 2 S„ 


co 2 kD 0 P G g I 3r 3 


d \ + 2^e + 


dr* 


SEg 

3r 


-k 





(23c) 


2 2 

where g' = dg/dr and g” = d g/dr . These equations must then be solved numerically. 
The electric and magnetic fields in the plasma are denoted as E r (r), Eg(r), and E z (r), 
and B (r), B„(r), and B (r). Applying boundary conditions gives the following expres- 
sion for the 9 - component of the electric field of an infinite- length current sheet at r = s: 


r 


[E^ ac (k)] _ =i^j* S K 1 (ks)^I 1 (ks) + 


Ki(ks) 


r 3E, 


I x (hp) I — — 1 - kI’ 1 (kp)(E 0 ) 


3r 


r=p 


r=P 


(24) 


r=s 2 
c 


'8E 


kK’ 1 (kp)(E 0 ) - K 1 (kp) ( 

r=p \ 3r 


9 1 


r =P 


with the additional constraint that 


A = 


-co 


kc 


Ki(kp)(E ) +K 0 (kp)/^)(B ) =0 

r=p \c / p r=p 


( 25 ) 


This additional constraint comes from the E z and Bg boundary conditions. Since E z (r) 
and Bg(r) are functions of Eg(r), equation (25) acts as a boundary condition for Eg(r) 
at r = p. 


Because of the particular numerical technique (see appendix C) used to solve for 
Eg(r), it is necessary to know the value of Eg(r) and its first three derivatives at r = 0. 
Because |Eg(r)J is zero, and Eg(r) is asymmetric about r = 0, Eg(r)] = 0 and 

[aXw/sr 2 ] = 0. From equation (24) it can be seen that only the ratio 



is required to calculate (EY ac ) . Since only the ratio 

v u 'r=s 

is needed and the differential equation for Eg is linear, only 
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the ratio |[9E 0 (r)/3r]y / (3 3 E 0 /3r 3 )| needs to be specified. Then, choosing 

f3E fl (r)/3rl =1 leaves only f3 3 E fl (r)/3r 3 to be determined by the following 
L ° J r=0 L J r-0 

iteration process. For r near zero, the terms in equation (23a) that involve the density- 
distribution function are small, and the solutions behave like Bessel functions. Hence, 

3 3 E fl (r)/3r 3 is first computed for the case of a constant-density plasma n = n . 

L 0 J r=0 

This value is then used as a first guess for the variable- density value. The values of 
E„(r), E (r), and E (r) at p are then computed, and the test of equation (25) is applied. 
If A * 0, a Newton- Raphson formula is used to find a new value for |3 Eg (r)/3r j and 

A is again computed. This iteration continues until the change in [3 3 E fl (r)/3r 3 ] from 

r 3 3T r =0 

one iteration to the next is must less than |3 Eg(r)/3r J . 

Once the values of [E 0 (r)J and [3E 0 (r)/3r] have been determined under the 

constraint A = 0, equation (24) can be used in equations (16). Now the same procedure 
used in deriving equations (17a) to (19c) may be applied to the variable- density case. The 
new expression for is 


R n = < 




! [K l( ks)]' 


Ii (kp) 



kI' 1 (kp)(E 0 ) 


r =P. 


r=P 


dk 


(26) 


k=k n 

r 3=° 


where 


r 3 (k) = kK' 1 (kp)(E 0 ) 


- K x (kp) 



(27) 


RESULTS AND DISCUSSION 

From the previous section, the independent variables in an ion cyclotron wave ex- 
periment included sheet-current radius, wavelength, and overall length; plasma radius, 
density, and radial density profile; steady magnetic field; and generator frequency. 
Also important are other effects, such as axial variations of the magnetic field, ion and 
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electron temperature, and damping phenomena. However, these effects have not been 
included because they entail more complex theory, and their effects can be minimized 
through proper design of the experiment. Conditions for which temperature and damping 
effects can be neglected are discussed briefly in this section. The intent here is to pre- 
sent results for a theoretical model of ion cyclotron wave generation that can be reason- 
ably approached in the laboratory. 


Current-Sheet Representation 

Figure 2(a) shows a Stix-type (ref. 1) radiofrequency coil that may be used to pro- 

i(k 0 z-u>t) 

duce a current sheet that varies approximately as j*e , where kg = 2ir/X. If N 

is the number of turns in each of the four sections of the coil, I is the peak radiofre- 
quency current flowing through the coil, and A q /M is the length of each section, the 
equivalent rectangular current sheet can be represented as shown in figure 2(b). The 
rectangular current sheet can be represented by a Fourier series with coefficients given 

as |A_„| = 4NM I sin where m = 1, 3, 5, etc. Since power transfer is proportional 
m mvA^ M 

i 1 2 u 

to |A m | , the fundamental (m = 1) is usually sufficient for power calculations. Thus, 
the rectangular current sheet is represented as 



(a) Radiofrequency coil. 


Surface 
current, j* 







j* ■ NMI/A. 

f _ 0 











M ' ‘ ' * T 
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(b) Current-sheet representation. 

Figure 2. - Radiofrequency-coil configuration and equivalent current sheet. 




Substituting this last expression into equation (19b) yields the power transfer from the 
fundamental component of the Stix coil: 


P 


t 


27ra 2 


4NMI 

V^o 



Z R n S r ( V 


( 29 ) 


Equation (29) can be used to compute power transfer for either the constant- or variable- 
density case, and for any value of N, M, X Q , I, a, p, s, n g , u>, and Aside from 
some practical electric circuit problems that are discussed briefly in this section, it is 
obvious from equation (29) and from the expressions for (eqs. (20), (22), and (26)) 
how the values of N, M, a, and w will affect power transfer. Hence, typical values of 
N = M = 4, a = 45 centimeters, and co/2 n = 6.5 megahertz are used for all calculations. 
The coil is 2 wavelengths long. Although the effects of coil wavelength and radius are ex- 
amined in this section, for illustrative purposes most of the results are presented for 
s = 10 centimeters and X Q = 45 centimeters. The aforementioned values of the parame- 
ters N, M, a, and u> are used in equation (29) to obtain 



Figure 3. - Radiofrequency-coil- 
radius correction factor as 
function of wave number for 
n th plasma mode. 


t = 1. 17X10 15 / 


watts 


lamperesy 


Z R nW 


n 


statohms 


(30) 


The effects of each independent variable on power transfer 
are discussed in the following paragraphs. 


Coil Radius 

The effect of coil radius s on P t is given entirely 
by the term 


F(k n s) = (k n s) 2 K 1 (k n s) 2 


which appears in (eqs. (20), (22), and (26)). This fac- 
tor F(k n s), is presented in figure 3 as a function of k n 
with s as a parameter. The term R^ (for s = 10. 0 centi- 


19 



Relative power transfer, P* statohm 


meters) may be corrected to any value of the coil radius by multiplying (for s = 

10. 0 cm) by the ratio F(k n s)/F(10k n ). However, approximate corrections to P^ can be 
made simply by multiplying P t (for s = 10. 0) by F(k o s)/F(10k Q ). It may be seen from 
F(k n s) that, for a given coil wavelength, it is desirable to make s as small as possible 
but subject to the restriction s > p. 


Frequency 

The function R^ is also directly proportional to u> (see eqs. (20), (22), and (26)). 
However, to increase the power transfer to the plasma merely by increasing the fre- 




Figure 4. - Relative power transfer as function of frequency ratio. 
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quency requires redesign of the coil system, since the radiofrequency coil impedance in- 
creases and causes high electrostatic voltages. Experimentally, the coil is operated so 
close to its natural frequency that further increases in co would introduce severe 
impedance- matching problems. In addition, larger, steady B-fields are required since 
peak power occurs at a given value of £2 = u)/S2.. 


Solid Cylinder of Constant-Density Plasma 


Figure 4 presents plots of relative power 


p* = 


o 2. *2 

2na. ] 


= Z R nW 


n 


as a function of f2 for all combinations of the following plasma radii and densities: 
p = 2. 5, 5.0, 7. 5, and 10.0 centimeters, and n g = 5x10^, 5x10^, and 5x10^ per cubic 
centimeter. An interpretation of these results and a discussion of significant trends fol- 
low. All the curves of P* as a function of Q, exhibit a peak (resonance). The value 
of £2 at the peak, as well as at the half power points, depends primarily on electron den- 
sity. For higher densities (n Q = 5x10 cm ), 12 at the peak is lower, and vice versa. 
Smaller dependencies on plasma radius and on current- sheet wavelength are discussed in 
more detail later. 

In general, the shape factor S r ( 7 ] n ) for the n = 1 mode controls the shape of these 
curves. The shape factor is given as a function of ri n in figure 5. 



Figure 5. - Shape factor S r (7j n ) as function of 7? n . 
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Plasma density. - An examination of figure 4 reveals that the optimum density for 

12 

power transfer to the plasma wave is about 5x10 per cubic centimeter or slightly less. 

This optimum density does not seem to depend strongly on plasma radius. In the section 

Variable Density and Annular Plasmas, it will be shown that power transfer is substan- 

tially lower than the maximum for n < 10 per cubic centimeter or n > 10 per 

cubic centimeter. 

Current- sheet wavelength . - The 4 5- centimeter wavelength chosen for most of the 
calculations does not necessarily represent the optimum wavelength for efficient power 
transfer. As a result, some calculations were made for current sheets 2 wavelengths 
long with wavelengths of 25, 35, and 55 centimeters. If it is assumed that I, N, and M 
are held constant, maximum power occurs when P* is a maximum. Tables I and II sum 
marize the results of these calculations by presenting P* and the corresponding 
value of S2 for all coil wavelengths, plasma radii, and densities for which constant- 
density calculations were made. 

Several features of table I should be explained, the first of which is that the optimum 
wavelength for any given combination of plasma radius and density depends on the actual 

TABLE I. - DEPENDENCE OF MAXIMUM RELATIVE POWER 
TRANSFER ON PLASMA RADIUS, ELECTRON DENSITY, 

AND CURRENT-SHEET WAVELENGTH 
[Surface current-sheet radius, 10.0 cm. ] 


Plasma 

column 

outer 

radius, 

P, 

cm 

Electron 

density, 

n e’ 

cm -3 

Wavelength of surface 

-current sheet, X Q , cm 

25 

35 

45 

55 

Maximum relative power transfer, P* „ . statohm 

9 max’ 

2.5 

5xl0 U 

5xl0 12 

5xl0 13 

0.87X10" 16 

1.97 

1.35 

0.97X10" 16 

1.98 

1.03 

0.84X10" 16 

1.50 

.612 

0.65X10" 16 

1.05 

.385 

5.0 

5X10 11 

3.08X10' 15 

3.24X10' 15 

2.64X10 -15 

2. 10xl0' 15 


5xl0 12 

5.03 

4.70 

3.50 

2.44 


5xl0 13 

3.08 

2.16 

1.24 

.726 

7.5 

5X10 11 

2.58X10" 14 

2.44X10" 14 

1.95X10 -14 

1.48X10' 14 


5xl0 12 

3.57 

3.09 

2.21 

1.49 


5X10 13 

2.41 

1.45 

.79 

.45 

10.0 

5X10 11 

1.32X10" 13 

1.08X10 -13 

0.80X10" 13 

0. 59xl0“ 13 


5xl0 12 

1.67 

1.27 

.85 

.56 


5xl0 13 

1.17 

.68 

.32 

.17 
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I 


sheet, * 0 , cm 


45 

55 

•atio, £2 = 

oj/£2. 

0.950 

0.930 

.80 

.75 

.43 

.36 

0.967 

0.953 

.81 

.76 

.43 

.37 

0.969 

0.958 

t-H 

00 

.76 

.43 

.37 

0.970 

0.958 

.81 

.76 

.43 

.37 


p* 

max 


table n. - ratio OF wave frequency TO values of p and n g . The optimum wave- 
length is shorter for high densities and large 

ION CYCLOTRON FREQUENCY AT MAXIMUM 

plasma radii. The optimum value of P^^ 
relative power transfer as function as a function of X Q is not always evident 

of plasma radius, electron density, from the table, since finding this optimum 

would require calculations for X„ less than 
AND CURRENT-SHEET WAVELENGTH O 

25 centimeters. For such short wavelengths 

Plasma Electron Wavelength of surface -current and for a given density, the value of Q, at 

column density, sheet, a , cm p* is closer to 1 so that it becomes in- 

outer n — IX13X 

e’ 25 35 45 55 creasingly difficult to satisfy the conditions 

FclCllUS j p “O 

D for which cold- plasma theory applies. How- 

p -' Frequency ratio, 12 = w/12. 

cm 1 ever, in determining the proper wavelength 

2.5 5X10 11 0.980 0.969 I 0.950 1 0.930 for the current sheet, other factors must be 

5xio 12 .91 .86 .80 .75 taken into account. One of these is ohmic 

| 5xio .62 .50 .43 . 36 losses in the current sheet (coil). For ex- 

5.0 5X10 11 0.987 0.978 0.967 0.953 ample, increasing the wavelength, and there- 

5xio 12 .91 .86 .81 .76 by increasing ohmic losses to get an in- 

13 

5xio .62 .51 .43 • 37 crease in power transfer for a given current, 

7.5 5X10 11 0.988 0.980 0.969 0.958 makes little sense if the change in ohmic 

5xio 12 .91 .86 .81 .76 loss exceeds the increase in power transfer. 

13 

5xio .62 .51 • 43 .37 Second, ion cyclotron wave experiments 

10.0 5X10 11 0.988 0.980 0.970 0.958 have been plagued by an undesired coupling 

5xio 12 .91 .86 .81 .76 of power to the plasma. This coupling was 

5xio .62 - 51 - 43 - 317 a result of high E„ fields in the plasma, 

which are produced by the voltages across 

the coil (ref. 3). This undesired coupling has been greatly reduced by the use of a prop- 
erly grounded electrostatic Faraday shield between the current sheet and the plasma. 

Since the voltage across the current sheet for a given current is determined primarily by 
the inductance of the current sheet, it may not always be desriable to change the wave- 
length (if it should increase inductance) to get increased power transfer. Furthermore, 
higher coil voltages can cause arcing problems between the coil and the ground points. 

Plasma radius . - Figure 6(a) shows the maximum value of P* as a function of 
plasma radius for several densities (X Q = 45 cm, s = 10. 0 cm). This power transfer 

P* increases approximately as the 4. 5 power of the plasma radius. Examination of 
m ax 

table I shows that this radius effect changes slightly with X . Figure 6(b) shows P* 

U IHcUv 

as a function of plasma radius for the case where s equals p (X Q = 45 cm). In this case 

P* , increases as the 3. 5 power of the plasma radius. These curves reveal that the 
max 

plasma radius has an extremely large effect on power transfer. 

Dispersion relation for finite electron mass . - For experimentalists, it is useful to 
have a simple expression relating the wavelength of the plasma wave to J2. Using basi- 
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Plasma 

Electron 

column 

density, 

outer 

V 

radius, 

cm -3 

cm 

P, 


cm 


2.5 

5X10 11 

5X10 12 

5xl0 13 

5.0 

5X10 11 

5xl0 12 

5X10 13 

7.5 

5X10 11 

5xl0 12 

5xl0 13 

10.0 

5X10 11 

5xl0 12 

5xl0 13 





cally the same model for the plasma as described in this report, but making the additional 
assumption that the mass of the electron was zero, Stix (ref. 1) derived the following dis- 
persion law that holds, approximately, for all plasma modes: 



(31) 


1 O 

For densities less than 10 per cubic centimeter and for plasma radii less than 5. 0 cen- 

th 

timeters, equation (10c) may be simplified. Let v j for the n 11 mode be denoted as 
If B 2 » 4C, then 


In 



o> 


and for m g * 0, 




A 2 / m 9 , 


(32) 


4tt 2 \ m i 


When boundary conditions were applied to the plasma model, the computer solutions 
showed that J 1 (^ 1 p) « 0 when r = 0. If the arguments of Jj.(i'jp) are denoted as X n 

9 9 

when T = 0, = (X n /p) . Thus, the relation obtained for the condition T = 0 is 



where X n = 3.83, 7.00, 10.15, 13.3, 16.5, etc., for n = 1, 2, 3, 4, 5, etc. 

The finite- electron- mass dispersion equation (eq. (33)) differs from the zero- 

electron-mass dispersion equation (eq. (31)) for small values of p, but the two equations 

2 2 2 

are identical in the limit of large plasma radius or high density (i.e., U*/c* » (X R /p) ). 
Equations (31) and (33) are plotted in figure 7. 

As noted before, the values of v^ n and A R (k n ) when T = 0 are the values obtained 
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:Zero electron massL 
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(a) Plasma wavelength, 20 centimeters. 
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(b) Plasma wavelength, 45 centimeters. 


Figure 7. - Frequency ratio as function of electron density for the natural mode with n = 1. 




ron mass 



Figure 7. - Concluded. 

for the natural modes of a plasma column surrounded by a vacuum. 

The computer calculations of Pj. showed that, for some values of p, n g , and ft, 
only the n = 1 mode makes a significant contribution to power transfer. However, for 
other values of p, n g , and ft, and especially for large values of p, several modes make 
a contribution. When several modes are excited in the plasma at one time, it may be dif- 
ficult for the experimentalist to make fundamental wavelength measurements to verify the 
existence of ion cyclotron waves. 


Variable Density and Annular Plasmas 

The solution for the variable-density cases required a greatly increased amount of 
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Radial density distribution function, g(r) 



Figure 8. - Radial density distribution as function of 
radius ratio. 


computer time; consequently, a more limited 
range of plasma parameters was studied. 
However, enough calculations were made to 
show the effects of radial density distributions. 
Figure 8 shows two density- distribution func- 
tions that were chosen to illustrate the effects 
of variable density: the bell-shaped and the 
parabolic. In addition annular plasmas of con- 
stant density whose inner radius, £ , was one- 
half the outer radius, p, were studied. Distri- 
butions similar to all three of these have been 
observed in ion cyclotron wave experiments to 
date. 

Table III gives P^. at maximum for a 
plasma radius p of 5. 0 centimeters, a 
current- sheet radius s of 10. 0 centimeters, 
and a current I of 92 amperes (I = 65 A) 

i lllo 

for several peak density values. The results 
are for the three density distributions and for 
constant density. Results for constant density 
from reference 2, where m g was assumed 
to be zero, are also presented in table HI. 

The effects of including electron inertia 


TABLE m. - POWER TRANSFER AT MAXIMUM AS FUNCTION OF 

PEAK DENSITY FOR SEVERAL DENSITY DISTRIBUTIONS 

[Plasma radius, 5. 0 cm; current-sheet radius, 10. 0 cm; current- 
sheet wavelength, 45 cm; current, 92 A.] 


Peak density, 

n o’ 

cm' 3 

Density distribution 

Constant, 

in = 0 
e 

Constant, 

m e^° 

Bell 
shaped, 
m e ^ 0 

Parabolic, 

m e^° 

Annular, 

m e t °> 

f =2.5 cm 



Power transfer, 

P t , kW 


5X10 10 

58.4 

8.2 

— 





5X10 11 

57.8 

26.6 

16.8 

23.2 

21.8 

5xl0 12 

42.4 

34.6 

33.6 

34.6 

31.6 

5X10 13 

13.0 

12.4 

23.8 

11.4 

17.8 
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3 -3 

terms may be seen by examining table m at high densities (n g =5x10 cm ) where the 
dispersion relations (eqs. (31) and (33)) are essentially identical. Here zero electron 
mass and finite mass give essentially the same value for P. for constant density. 

Lj lllciX 

At lower densities, however, the zero-mass results are too high (by an order of magni- 
tude or more for n g < 10 cm ). 

For the coil configuration of figure 2(a) (p. 18), there are two current branches. If 
the current I in each branch is 92 amperes, as was assumed for the calculations in 
table HI, the total current is 184 amperes. Swett (ref. 3) measured the total resistance 
of a coil of this type and found it to be about 0. 3 ohm. Thus, the ohmic loss for 92 am- 
peres is approximately 5 kilowatts. This analysis would then predict that power -transfer 
efficiencies for the best cases (m g * 0) in table HI could be as high as 87 percent. 

In table HI, the use of peak density as a parameter for comparing the different radial 
density distributions was arbitrary. A more significant parameter might be the volume 
average density given by 


n 


e 


2vn /*P 
/ g(r)r dr 

7Tp 2 


where n Q is the maximum density in the plasma. For the bell- shaped distribution, 
n g = 0. 312 n Q ; for the parabolic distribution, h g = 0. 653 n Q ; and for the annular distribu- 
tion, n_ = 0. 750 n . 

’ e o 

In figure 9 P* is shown as a function of this volume average density for the 

rricix 

constant- density case and for the three variable-density distributions. The points for the 
variable- density plasma are in good agreement with the curve for constant density when 
plotted in this manner. As expected from the strong dependence on plasma radius, 



Figure 9. - Relative power transfer as function of volume average electron density. 
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power transfer to the hollow plasma is almost as great as that to a constant-density 
plasma with the same outer radius and the same peak density. Apparently, the plasma 
nearest to the coil carries most of the wave energy. 

The importance of the low-density plasma in the tail of the bell- shaped distribution 
is seen by comparing P, for this distribution with that for the parabolic distribu- 

iii ax 

tion (in table m) for the same peak density. When the peak density is too high, 

13 -3 

(n =5x10 cm ) the bell shape is better; but when the peak density is too low 
° 11-3 

(n Q = 5X10 cm ), the parabolic shape is better. When the peak density is about opti- 
mum with respect to relative power transfer, P. is the same for the two distribu- 

t j max 

tions. The P. for a constant density plasma whose radius is 2. 5 centimeters (in- 

t y max 

stead of the 5. 0 cm assumed in fig. 9) is less than that of the bell- shaped distribution at 
all peak densities. Thus, the tail of the bell- shaped distribution is important at all peak 
densities, but more so at high densities. 

The shapes of the curves of P^. as a function of ft were the same for variable den- 
sity as those for constant density if the peaks occurred at the same value of ft in both 
cases. The widths of the curves are only slightly dependent on the number of modes for 
which there is significant power input; thus, the widths of the absorption maximums de- 
pend primarily on the value of ft at which peak power transfer occurs. 


Different Plasma Modes 

Each term in the sum in equation (19b) is regarded as the amount of power put into 

one of the natural modes of the plasma wave. If the plasma radius is small (<5. 0 cm) 

12 -3 

and the density low (<10 cm ), the values of k for the natural modes are widely 
spaced; so that if kj = kp, the majority of the power goes into the first mode because the 
shape factor S r (f? n ) (fig- 5, p. 21) is small for all k n 's with n > 1. When the plasma 
radius is large and the density is high, the k n 's are closely grouped, and there can be 
significant contributions to power for several modes. 

It might be expected that, at small plasma radius and low density (where k^s are 
widely spaced), a curve of P^ as a function of ft might show a maximum each time the 
k for one of the modes equaled k Q . This would be true if R n were the same order of 
magnitude for several modes. However, the calculations showed that R n increases or 
decreases with n in such a way that only a single significant peak results. This single 
significant peak was found for all constant-density calculations as well as for the variable- 
density calculations. Thus, if this is the only manner in which power is transferred to 
the plasma, an experimentalist should see only a single resonance peak in P^ as he 
varies the magnetic field. Of interest with respect to the results is that for large plasma 
radii and high density, R n is a complicated function of ft with many maximums and min- 
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imums. Because of the complexity in R n , the modes that absorb most of the power are 
continually changing as the magnetic field is varied. For some instances, it was neces- 
sary to calculate contributions from as many as 20 modes to determine the total power 
absorbed by the plasma. When the power is distributed among many modes, the struc- 
ture of the radial variation of the electric and magnetic fields can be quite complicated. 


Effect of Temperature and Collisions 

Since one of the reasons for the interest in the production of ion cyclotron waves is 
their potential for heating the ions of a plasma, a cold plasma model may seem inappro- 
priate. Even when the heating is supposedly localized at a magnetic beach removed from 
the wave-driving portion, some wave damping may occur throughout the plasma. The 
applicability of cold plasma results to a hot, collisionless plasma must, therefore, be 
examined. Stix showed that the dispersion relation for a hot, collisionless plasma re- 
duces to that of a cold plasma when 


m -(co - fi.) 2 

— — » 1 

2/cTdt 2 

For a plasma with density n e = 5x10^ per cubic centimeter, peak power transfer 
occurs at £2 ~ 0. 95 and 


rm(c o - £2^ mi2? (1 - £2)^ 

2/cTVk 2 2/<T i k 2 


100 

T^(electron volts) 


Thus, for T. < 10 electron volts, the cold- plasma dispersion relation is good at low den- 
1 11 

sities, n e = 5x10 per cubic centimeter. At higher densities, 1 - £2 is larger at reso- 
nance, and the cold-plasma dispersion relation may be used at even higher temperatures. 

Under certain conditions, a low- density plasma with a neutral background gas or a 
high-density plasma may become collision dominated. The following expressions for ion- 
electron and ion- neutral momentum transfer collision frequencies are given in refer- 
ences 4 and 5, respectively: 


^ie 


%4 
m i 3 



n e e 


( kT e ) 


3/2 


In 


3 (kT) 3 / 2 


2V ? t 


e 3 n y 2 

e 
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^in = 1x105 P’ 

where p' is the neutral gas pressure in microns. For an electron temperature of 1 elec- 
12 

tron volt, n g = 5x10 per cubic centimeter, and p' = 2 microns, the ion- electron colli- 
sion frequency is 


v. = 8x10^ sec 

the ion-neutral collision frequency is 


v in = 2xl0 5 sec" 1 


and the ion collision frequency is 

u i = v ie + ^in = 0x105 se c -1 

4 

For an electron temperature of 10 electron volts, is even less than 8x10 per second. 

It is probably sufficient to say that the plasma is collisionless when v. « co/2n so 
that each ion makes many orbits between collisions. The value chosen for u>/2n for the 
calculations in this report was 6. 5x10° per second. Since many other ion cyclotron wave 
experiments in the future may be run at even higher frequencies, the plasma does not be- 
come collision dominated unless the density is high and electron- ion collisions predomi- 

12 

nate. In fact, at the densities of 5x10 per cubic centimeter obtained in this study for 
optimum power transfer (see tables I and n), the plasma may be considered collisionless. 


Wave Reflections 

The calculations in this report were made by assuming the presence of an infinitely 
long plasma column immersed in a uniform magnetic field. Actual laboratory plasmas, 
however, do not satisfy this assumption, in that, at some distance from the wave- 
generating current sheets, there may be magnetic mirrors, magnetic beaches (regions 
where O « 1), conducting or nonconducting end walls on the plasma test section, and/or 
other nonuniformities. It is more than likely that some of these nonuniformities will be 
the cause of wave reflections from the region of nonuniformity back in the direction of the 
generating structure. If the ion cyclotron wave is sufficiently damped in these regions, 
the reflected wave is small and its effect on power transfer is negligible. However, a 
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large reflected wave will pass back under the wave- generating current sheet and may 
cause constructive or destructive interference there, thus altering the power transfer. 

To compare this theory with experiment, reflections should be eliminated. It is expected 
that the most serious wave reflections come from magnetic beach regions, where for low 
ion temperatures the index of refraction k n c/w changes most rapidly for only small 
changes in magnetic field. To study these reflections, it is necessary to take into ac- 
count the effects of nonuniform magnetic fields; near S2 = 1, hot plasma theory must be 
used. 


CONCLUSIONS 

Calculations of the power transferred to ion cyclotron waves in both homogeneous 
and inhomogeneous plasmas by a Stix coil have been made. These calculations show 
that, for a fixed current- sheet geometry and a constant- density plasma, the power trans- 
fer is roughly proportional to the 4. 5 power of the plasma radius. For equal plasma 
current- sheet surface and column outer radii, the power varies as the 3. 5 power of 
plasma radius (or coil radius). If the total current in the current sheet is assumed con- 
stant, there is an optimum wavelength that depends on plasma density and radius. The 

12 

optimum volume average electron density for power transfer is about 5x10 per cubic 
centimeter. 

12 

For electron densities less than 10 per cubic centimeter and for plasma radii less 
than 5. 0 centimeters, most of the power is absorbed by the lowest natural radial mode 
of the plasma. For electron densities greater than 10 13 per cubic centimeter and for 
plasma radii greater than 7. 5 centimeter, the power is distributed between several nat- 
ural modes with the lowest mode not necessarily absorbing the most power. For all 
cases in which calculations were made, a plot of power transfer as a function of magnetic 
field gave only one significant resonance peak, regardless of how many modes were ab- 
sorbing a significant amount of power. 

Comparison of the calculations of this report with those of a previous report (NASA 

TN D-3361), in which a zero- electron- mass dispersion law was used, shows that the 

zero- electron-mass assumption is valid for a particular mode when the plasma electron 

12 

density is high. For an electron density less than 5x10 per cubic centimeter, the 
zero-electron-mass calculations give significantly greater power transfer than that cal- 
culated in this report. 

Calculations for four different radial density profiles show that power transfer to in- 
homogeneous plasma columns is high when the volume average density of the plasma is 
near that density for which power transfer is optimum for constant-density plasmas. 
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Finally, the applicability of the results presented in this report depends on the degree to 
which a laboratory plasma is represented by the model used herein. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, December 2, 1966, 

129-01-05-09-22. 
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APPENDIX A 


SYMBOLS 


B 

B 

B(r) 

B„ 


B plasma(vac) 

B r 

^.plasmafvac) 

B z 

R plasma(vac) 

a 9 


c 

D 


one-half current-sheet length 

k 2 (P + S) _ w? p + (D g - S 2 ) 

S c 2 c 2 S 

magnetic field vector (time-varying part) 

radial dependence of B 

magnitude of B Q 

superimposed static magnetic field parallel to z-axis 


radial component of B in plasma (vacuum) 
axial component of B in plasma (vacuum) 
azimuthal component of B in plasma (vacuum) 


P 

S 


k 4 - 2 — k 2 S + — (S 2 - D 2 ) 

2 4 

c c 


velocity of light 
D' 


D' 

D o 

E 

T-,plasma(vac) 

"r 

^plasmafyac) 

pplasma^ac) 



e 


F(k n s) 

g(r) 

I 


eq. (5g) 

D at peak density 
electric field vector 

radial component of E in plasma (vacuum) 

axial component of E in plasma (vacuum) 

azimuthal component of E in plasma (vacuum) 
azimuthal component of finite- coil electric field 

electronic charge 

radiofrequency-coil-radius correction factor 

radial density distribution function 

peak radiofrequency current in current sheet 


Iq modified Bessel function of first kind (zero order) 

Ij modified Bessel function of first kind (first order) 

Jq Bessel function of first kind (zero order) 

J j Bessel function of first kind (first order) 

j* magnitude of j* 
j* surface current density vector 

j(r) plasma current density vector 

K plasma dielectric tensor 

K' K- _1 

Kq modified Bessel function of second kind (zero order) 
modified Bessel function of second kind (first order) 
k plasma wave number 

J.U 

k n wave number for n m plasma mode 
k Q current- sheet wave number 
L L’ + 1 

L’ eq. (5d) 

& linear differential operator, — ( — — r) 

0r l r dr J 

S. plasma column inner radius 

M fraction of wavelength covered by one section of Stix coil (see fig. 2) 
M tensor, T(k) = det M for annular plasma column 

Mjj elements of M 

m g electron mass 

m i hydrogen ion mass (proton mass) 

N number of turns in one section of Stix coil 

n(r) electron (ion) density 

n g electron density 

n g volume average electron density 

n. ion density 

n Q peak electron density 
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p P’ + 1 

P Q P at peak density 

Pj. power transfer (see eq. (19a) 

P' eq. (5e) 

2 *2 

P* relative power transfer, P^./27ra j 

P* „ v maximum value of P* on plot of P* against J2 

IIlaA 

p plasma column outer radius 

p' background neutral gas pressure 

R R’ + 1 


R’ 

R n 

R n 

r 


S 

S’ 


W 



v e (r) 

v^r) 



z 


a 




eq. (5c) 
eq. (17b) 
eq. (19e) 

radial cylindrical coordinate 
S' + 1 


eq. (5f) 

S at peak density 

shape factor (see eq. (19c)) 

surface current- sheet radius 

electron temperature 

ion temperature 

time 

electron vector velocity 
ion vector velocity 
Bessel function of second kind 
Bessel function of second kind 



kD 


kD 


(zero order) 
(first order) 
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g g(k) denominator of eqs. (20), (22), or (26) 

! KjOcp) 

1 kp K Q (kp) 


A 


^n 

9 

A 

9 


K 


*n 


v ie 


in 


v 


2 

1 


v 


2 

2 


n 


e 


n i 

a 

w 


1 


boundary condition (see eq. (25)) 

2(k Q ' k n> a ( eq ‘ ^ 19d ^ 
azimuthal cylindrical coordinate 
azimuthal unit vector 
Boltzmann constant 

it. 

wavelength of n cn plasma mode 
wavelength of surface current sheet 
ion collision frequency 
ion- electron collision frequency 
ion- neutral collision frequency 

-B + Vb 2 - 4C 
2 


-B 


Vb 2 - 4 C 


electron plasma frequency 


, (4 7 rn e e 2 / 4i e ) 


1/2 


, ^n. 


e 2 /m. 


\l/2 


i ) 


ion plasma frequency 

w/Slj 

electron cyclotron frequency, eB Q /m e c 
ion cyclotron frequency, eB Q /m.c 
angular frequency of plasma wave 
unit matrix 
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APPENDIX B 


FUNCTIONS T(k) AND MATRIX M 

The functions Tj(k) (eq. (20)) and the matrix M used for computing r 2 (k) are given 
as follows: 

For the solid cylinder (£ = 0): 

ryk) = afpjjCi'jp) + yJ'jpJ^i'jP)] [yy 2 pJ 0 (y 2 pJ 0 (^2P) + 

- /sjpj^p) + yt' 2 pJ 0 (t’ 2 p)j [yi' 1 pJ 0 (i' 1 p) + 

For the annular cylinder (£ ^ 0): 

M = (M ij ) i,j = l,2,3,4 

r 2 = det M 
M 11 = “ J l^i £ ) “ 

M 12 = P J &2^ - 

M 13 = a - 

M 14 = 0 Yfal) - 


Vl J 0 ( Vl H)I l( M) 

Pkl 0 (k£) 

PkI 0 (M) 

Pkl 0 (k£) 

Pkl 0 (k£) 


^21 - 


^(^^(kt) 

kl 0 (k£) 
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M 22 = 


VoM 1 ^ 

kI 0 (M) 


M 23 " Y lM 


^oM^) 

kl 0 (kf) 


M 24 = Y 1<V> " 


kIo(k£) 


M 31 = a J^p) + 


*' 1 J 0 ( l ' 1 P) K l< k P> 

PkKjjOcp) 


M 32 = ^ J 1^2 p) + 


^ 2 J 0 ^ 2 p ^ K 1 ^ 

Pkl^Ocp) 


M 33 = a Y l^lP) + 


^ Y 0 (^p) K i(kp) 

PkK^kp) 


M 34 = ^ Y 1 ^ 2 p ) + 


PkKQ(kp) 


M 41 = + 


^l J 0 (^i p ) K l(k p ) 

kK()(k p ) 


M 42 = J 1^2 p) + 


kK Q (kp) 


M 43 = Yj^jp) + 


,l Y 0 (.i p ) K i(kp) 

kKQ(kp) 


M 44 = YjC^p) + 


^ 2 Y 0 ^ 2 p ^ K 1 ^ 

kKfjCkp) 
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o 

In the preceding expressions for and M.j, 

When or is negative, ^ is replaced by 

Ij(x), Yq(x) by -Kg(x) and Yj(x) by K^x), where 
or 4. However, a and /3 are not changed. 


and were assumed to be positive. 

v x > by j i( x ) by 

x = y| 2 1 * y> and y is either p 
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APPENDIX C 


SOLVING FOR NATURAL MODES 

For the constant- density plasma where the electric fields of an infinitely long plasma 
column are known analytic functions of k, the condition r = 0 for natural modes results 
in one equation with k as the only variable. Then, k is varied until this equation is sat- 
isfied. 

For the inhomogeneous plasma, the boundary conditions can be combined so that 
there are two conditions, T = 0 and A = 0, for natural modes, and also two unknowns, k 
and |a 3 E 0 (r)/ar 3 J ^ = E"\ Thus, 


ar ar 

r=r + — - dk + -dE ,M (Cl) 

° ak aE"' 


where 


3A 3A 

A = A + 5 dk + 2. dE"’ 

3k aE**’ 


r ' r < k b’ E b") r o= r < k a > E a'> 

A = A < k b- E b"> A D =A(k a ,E k ") 

dE f " = Ej," - E'" dk = - k a 


(C2) 


The values of r Q and A q are the values of T and A for a set of trial values of k 
and E ,M , to be designated k and E'". Setting r=0 and A = 0 in equations (Cl) 
and (C2) yields two equations and two unknowns, dk and dE M * that may be determined 
by Cramer's rule: 
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dk = 


0 
(-( 

1 

9r o 

0E ,M 


ar o 

0k 

1 

*1 

o 

0 
« 

1 

9A o 

0E ,M 

dE™ = 

9A o 

0k 

— 

0 
< 

1 

9r 0 

ar o 


9r o 

CD 

O 

9k 

0E ,M 


0k 

0E ,M 

9A o 

o 

<1 

CO 


CD 

> 

o 

9A o 

9k 

0E ,M 


0k 

0E ,M 


The value of k is then changed by dk and the value of E' M is changed by dE'"; and 

a 

A q and r Q are computed again. This Newton-Raphson procedure is repeated until A Q 
and r Q converge sufficiently close to zero. 

Two problems arise when the previous procedure is used: first, because of the na- 
ture of r and A, it is possible that there will be no convergence at all; second, there 
may be convergence, but to the wrong mode. These problems can be overcome by plot- 
ting r as a function of k, under the single constraint A = 0 for a single set of the 
values, 0 and 11^. The T = 0 points then give the values of k for the natural mode. 
Since this latter technique takes considerabe computer time, the following procedure was 
adopted: 

(1) For a given density distribution g(r), values for and IT^ were chosen and 
the procedure described in the previous paragraph was followed to find the values of k 
and E'” for the natural modes. 

(2) Trial values of k for other values of n„ and were obtained from the 

a a e a 

cold- plasma zero- electron-mass dispersion law and from the values computed in item (1): 


k 


2 

a 


k 


2 


« nj q - „») 

n e (‘ - n a) 


(3) These trial values of k & obtained in item (2) and E”’ from item (1) could then 
be used in the Newton-Raphson method to find natural modes. 
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